Supersymmetry at LHC and ILC* 



K. Desch 

University of Hamburg, 22761 Hamburg, Luruper Chaussee 149, Germany 
now at University of Freiburg, 79104 Freiburg, Hermann-Herder-Str.3, Germany 



The prospects for the discovery and exploration of low-energy Supersymmetry at future colliders, the Large Hadron 
Collider (LHC) and the future international linear electron positron collider (ILC) are summarized. The focus is on 
the experimental techniques that will be used to discover superpartners and to measure their properties. Special 
attention is given to the question how the results from both machines could influence each other, in particular when 
they have overlapping running time. 



1. INTRODUCTION 
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The search for SUSY and, should it be found, measurements of the supcrpartner properties are among the most 
' important motivations for future high energy particle colliders. A general introduction to TeV-scale Supersymmetry 
<N ! (SUSY) as one of the best motivated extensions of the Standard Model (SM) has been given elsewhere in these 
proceedings The Large Hadron Collider (LHC) currently under construction at CERN will go into operation 
t-H 1 in 2007 and has a huge potential for SUSY discovery as well as for first measurements of SUSY particle properties. 
■ The planned International Linear electron positron Collider (ILC) is an ideal tool for precision SUSY measurements. 
, Both machines together will be able to give important insight into the mechanism of SUSY breaking and may open 
' a window to GUT/Planck scale physics. 

In this article, in Section the prospects for inclusive SUSY discovery at the LHC arc summarized. Furthermore 
I/"} . studies for the exclusive reconstruction of superpartners and the measurement of their masses are explained. In 
| Section |3 the prospects for precision SUSY measurements at the ILC are shown. The different techniques for mass 
measurements, measurements of polarized cross-sections and quantum numbers of the superpartners are detailed 
^p" 1 - with some explicit examples. Finally, in Section 01 the interplay of the anticipated results of both LHC and and ILC 
Qh. in particular when analyzed simultaneously is shown. The extended Higgs-sector of supersymmetric models also is 
an integral part of future exploration of SUSY. The prospects for Higgs searches and precision measurements at LHC 
^ ' and ILC are summarized elsewhere in these proceedings 0. 

•rH , 
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c3 2. SUSY AT THE LHC 

2.1. Experimental environment 

The LHC will produce proton-proton collisions at a center-of-mass energy of 14 TeV initially at a luminosity of 
1 x 10 33 cm _2 s _1 (low-luminosity) and later at 1 x 10 34 cm _2 s _1 (high luminosity). The high center-of-mass energy 
makes this machine well suited for the direct discovery of new massive particles beyond the SM. The partonic 
luminosity is large enough to pair-produce colored new particles (like squarks and gluinos) up to masses of a few 
TeV at observable rates. The high luminosity and the large total pp cross-section impose strong requirements on 
the performance of the detectors and in particular on the trigger. The two multi-purpose detectors, ATLAS 0] 
and CMS 0, will be able to trigger efficiently on high-pT leptons, jets and are due to their hcrmcticity sensitive to 
missing transverse energy. Thus they are well suited for the typical experimental signatures of SUSY. The hadronic 
environment, however makes the exclusive reconstruction of final states as well as precision measurements challenging. 
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Figure 1: SUSY mass spectrum for the mSUGRA benchmark point SPSla 



This is due to generally huge QCD backgrounds, the presence of pile-up events at high luminosity, and the absence 
of a longitudinal beam constraint. In the past years several sophisticated analysis techniques for SUSY processes 
at the LHC have been developed which in spite of the difficult conditions go far beyond inclusive SUSY discovery 
and will allow for a significant set of first measurements of superpartner masses and some of their properties at least 
under favorable circumstances. 



2.2. Production processes 



At the LHC, the predominantly produced superpartners will be the colored gluinos and squarks. If R-parity is 
conserved, they will be pair-produced at large rates (typically O(10 pb) at masses around 1 TeV), comparable to the 
SM jet rates at the same values of Q 2 = Mg USY - Direct production of sleptons, charginos and neutralinos mainly 
proceeds via Drell-Yan production and t-channcl squark exchange at a much lower rate. However, the color-neutral 
superpartners often appear in the decay chains of squarks, if kincmatically allowed. 

If R-parity is conserved, the Lightest Supersymmetric Particle (LSP) is expected to be neutral, stable and only 
weakly interacting. It escapes detection and leads to the most distinctive SUSY signature: large missing transverse 
energy since all superpartner decay chains eventually end in the LSP. In mSUGRA and AMSB SUSY breaking 
models, the LSP is the lightest ncutralino, in GMSB it is the gravitino. 

In Fig. the superpartner spectrum for typical mSUGRA point, SPSla Q is shown. This particular benchmark 
scenario has been extensively studied both for the LHC and for the ILC. It provides a very rich phenomenology at 
both machines since the complete spectrum lies below 600 GeV. For the LHC also a larger set of post-LEP benchmark 
points corresponding to various SUSY breaking mechanisms and parameter sets has been studied. 

If R-parity is broken, the missing energy signature gets lost. Depending on the R-parity violating model multi-jet 
and/or multi-lepton signatures arise. They have also been studied 0] but will not be further discussed here. 



2.3. Inclusive discovery 



Due to the large production cross-sections, the SUSY particles can be inclusively observed over the SM background 
in the LHC data with very simple cuts. The generic signatures are large missing transverse energy (fix) and multiple 
hadronic jets and/or leptons. A typical example of this signature is the distribution of the so-called effective mass, 

4 

M cS = p T + J2PT,i, 

i=l 

i.e., the sum of the missing transverse energy and the transverse energy of the four hardest jets. Its distribution is 
shown in Fig. [21 together with the expected SM background for a mSUGRA model with squark masses of approxi- 
mately 700 GcV after requiring at least four high-P^ jets and significant missing transverse energy. The simulated 
data correspond to an integrated luminosity of 10 fb _1 . It can be seen that for large values of M e g SUSY events can 
be selected with negligible SM background. 
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Figure 2: M e g distribution for a typical mSUGRA point and SM backgrounds after cuts. 

In Fig. |3fleft) the 5tr discovery reach in the plane of the mSUGRA parameters m and m 1 / 2 is shown for inte- 
grated luminosities of 1,10,100,300 fb" 1 (red lines). The remaining parameters Aq and tan/3 are fixed to and 35, 
respectively, and sgn(/i) is chosen positive. The discovery reach depends only weakly on these parameters. Also 
shown arc lines of constant squark and gluino masses. For (1,10,300) fb _1 the mass reach for squarks and gluinos is 
approximately (1,2,2.5-3) TeV thus covering a very large part of the mSUGRA parameter space. The same applies 
qualitatively as well for a large part of the general MSSM with neutralino LSP. 

Under the assumption of GUT unification of the gaugino mass parameters Mi, M2, M3 and sfermion mass param- 
eters, the gluinos and squarks are heavier than the color-neutral supcrpartners due to renormalization group effects. 
Hence, the heavier charginos and neutralinos often appear in the decay chains of squarks. Their electro-weak decays 
often give rise to high-p^ leptons whose efficient detectability provides an additional inclusive SUSY signature in 
many cases. In Fig. fright) the reach of the different lepton signatures (1 lepton, 2 like-sign (SS), 2 opposite-sign 
(OS), 3 leptons) plus ,Et is shown for 10 fb _1 . 




Figure 3: Left: LHC reach in the mo — nii/2 plane for inclusive discovery of SUSY in the missing transverse energy signature 
for the CMS experiment for an mSUGRA scenario with tan j3 = 35. Right: LHC reach for various inclusive SUSY signatures 
involving 0-3 leptons for 10 fb . 



2.4. Measurement of SUSY particle masses 



While inclusive detection of SUSY processes is quite straight-forward at the LHC, the reconstruction of exclusive 
decay chains and the reconstruction of superpartner masses is quite involved due to various reasons: 1. the long 
decay chains of gluinos and squarks lead to signatures with many jets and leptons with huge combinatorics. 2. Due 
to the unknown longitudinal boost of the colliding partons of the initial state no kinematic constraints from the 
beam particles can be applied. 3. The event-by-event reconstruction of superpartner invariant masses is not possible 
due to the undetected LSP's. The reconstruction of superpartner masses has therefore to rely on the detection of 
kinematic endpoints in the invariant masses of the detectable final state partons (jets and leptons) as well as on some 
knowledge about the involved decay chains. 




Figure 4: Typical decay cascade of a left-chiral squark 



2.4.1 . The q L — ► q£ + t Xi decay chain 

A frequently occurring and rather well reconstructible decay chain of a left-chiral squark via the second-lightest 
neutralino x® an d a right-chiral slepton £r is shown in Fig. 0] The invariant mass distribution of the two opposite- 
sign same-flavor (OS-SF) lcptons has a characteristic triangular shape which exhibits a distinct kinematic endpoint 
('edge') which involves the unknown masses of the three involved SUSY particles: 

M^ e = V /(M|-A/|)(M?-A/| ) )/A^. 

In Fig.OJlcft) the (OS-SF) di-lcpton signal from x§ decay is shown is shown together with the background from other 
SUSY decays and the (negligible) SM background. The SUSY background results mainly from wrongly combined 
and thus uncorrelated lcptons from independent neutralino or chargino decays. It can be efficiently determined from 
the rate of OS-OF di-lcptons in the data and then subtracted, as shown in the right part of Fig. 03 If the sleptons 
are heavier than the three-body decay ~ * i + £~Xi dominates which is different in shape and has an endpoint 
at the neutralino mass difference, A/^o — Af^o. 

' Ao AT 




Figure 5: Left: opposite sign same flavor (OS-SF) di-lepton mass spectrum. The upper curve is from X2 decays, the next curve 
is SUSY background, the lowest curve is SM background. Right: same after subtraction of the OS-OF rate which reduces the 
background from final states with uncorrelated lepton flavour. 

While from the di-lepton endpoint alone no absolute superpartner masses can be extracted, further information can 
be obtained from the squark decay chain shown in Fig. 0] from various combinations of the lcptons with a jet 
In particular, the following additional mass relations for kinematic edges can be exploited: 
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Figure 6: Various mass spectra which exhibit kinematic endpoints in squark decay cascades. 



The labels "min" and "max" refer to the distribution constructed from the smaller and the larger of the two q£ masses. 
Furthermore "thres" refers to the threshold in the subset of the m q u distribution for which the angle between the 
two lepton momenta (in the slepton rest frame) exceeds tt/2, which corresponds to mf, se /V2 < mu < m^ se . The 
corresponding mass distributions are shown in Fig. El The position of the di-lepton edge can be measured to a 
statistical precision of better than 100 MeV and the edges involving jets can be measured to few GeV precision with 
100 fb _1 of data. In the case of the latter a systematic uncertainty from jet energy scale of approximately 1% has to 
be accounted for as well. Since the number of measurable edges in this scenario is larger than the number of involved 
superpartner masses, the absolute masses can be extracted from a simultaneous fit. The achievable precisions are 
listed in the first column of Table [n] 



2.4.2. Gluino and third generation squarks 

Knowing the masses of XiiX2^R' an d Ql, gluinos can be reconstructed by adding another quark jet to the q££ 
system. In particular when the decay proceeds through a b squark, the tagging of two b-jets reduces combinatorial 
background. Knowing the mass of the LSP (from the joint fit to the kinematic edges described above), the momentum 
of the X2 can be approximated by 

m(XiY 



m(££) 

if the \i carries negligible momentum in the x® res t frame, which in the SPSla scenario is the case for events close 
to kinematic endpoint of the di-lcpton mass spectrum. Knowing the X2 momentum and mass, the sbottom mass 
can be subsequently reconstructed as the 6^2 invariant mass and the gluino mass as the 66x2 invariant mass. The 
reconstructed mass distributions in CMS for sbottom and gluino are shown for 300 fb _1 are shown in Fig. (left, 
middle) . The mass difference between gluino and sbottom can be reconstructed without assumptions on the Xi mass 
(Fig.|7| right) i. 

The reconstruction of t squarks is more challenging. Initial studies are available 
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Figure 7: Reconstructed mass distribution for sbottom (left), gluino (middle), and for the gluino-sbottom mass difference 
(right) (from %). 



2.4.3. Mass relation method 

An alternative approach to mass reconstruction in the presence of long decay chains at the LHC, the mass relation 
method has recently been proposed (icj . If one considers e.g. the decay chain g — > bbi — > — ► £b\bi£i — > 
Xi^i^2^i^2, the five superpartner mass m^o, m^, to^o , nig can be calculated from the 5 4-momenta of the final 
state particles all of which except for p%o are measured. Thus, the set of superpartner masses compatible with a 
single observed event corresponds to a 4-dimcnsional hypersurface in the 5-dimensional mass space. Since the exact 
location of the hypersurface is different for each event, the ensemble of hyper-surfaces from all events will have 
intersections at the true values of the five unknown masses. This method has been applied in a simplified version to 
the reconstruction of sbottom and gluino mass under the assumption that the other masses are known, in which case 
the hypersurface reduces to a line in (mg,mj;) space, i.e. for each pair of events a solution for (mg,m^) is obtained 
up to a two-fold ambiguity. The advantages of the method are that it is based on exact kinematics without any 
approximation and that all events, not only the ones close to the kinematic endpoint, can be used. Furthermore, 
mass peaks are reconstructed rather than kinematic edges. 

2.4.4. Sleptons 

If sleptons are sufficiently light, they arc produced at a decent rate in pp collisions through Drcll-Yan pair pro- 
duction. For the SPSla benchmark, the masses of the left- and righthanded selectrons and smuons are 143 GeV and 



202 GeV, respectively. The total cross-section for smuon/selectron pair production is 91 fb. The signature are two 
opposite sign same flavor leptons, missing transverse energy and no jets. After subtraction of the opposite flavor 
background only a few events remain for 100 fb -1 . However, from a transverse mass estimator, the slepton mass can 
be estimated to a few GeV precision under favorable circumstances [llj . 

Staus may frequently occur in the decays of xH if the fx is lighter than X21 which is the case in SPSla and in 
generally in mSUGRA models with significant f mixing at large tan/3. Hadronic r decays can be tagged typically 
at an efficiency of 50% for a QCD jet rejection factor of 100 at low luminosity. Selecting di-r events with fix and 
large M e s and again subtracting the same sign contribution, the invariant mass of the rr decay products carries 
information about the f mass in its endpoint. A mass estimate with a couple of GeV precision seems feasible but 
further study is needed. 

2.4.5. Charginos 

Recently a method has been proposed to observe the lighter chargino, xti frequently appearing in the decay chain 
of a left squark, — > xfq via its decay x^ — * W^Xi — ► qq'xi- The methods uses di-lepton events where the two 
leptons arise from the decay of the initial squark decaying into X2 as described above. Assuming known masses for 
XiiXti^Ri an d <7l, the momentum of the Xi from the x° can be reconstructed up to a two- fold ambiguity. After 
identification of hadronic W bosons, the mass of the chargino can be reconstructed from the W momentum, the 
reconstructed opposite-side Xi an( i the total JEr up to another two-fold ambiguity. After background subtraction 
of events in the side-bands of the reconstructed W mass, a peak becomes visible in the reconstructed xt mass 
distribution. For a mSUGRA scenario with mo = 100 GeV, m 1 / 2 = 300 GeV, Ao = —300 GeV, tan/3 = 6, and 
sgn(^) = +, a 3a excess is achievable for 100 fb -1 and a mass estimation with approximately 10% error seems 
feasible Q . 

2.4.6. Heavy gauginos 

The left squarks predominantly decay via the gaugino-like neutralinos and charginos, i.e. usually X2 an d xf m 
mSUGRA models. While xi often is pure Higgsino, the heaviest neutralino xi an d the heavy chargino xf nave 
some gaugino admixture leading to a production of xi arL d xf m the q~L decays with a branching ratio of a few 
percent. Four different decay chains via left and right sleptons (for the neutralino) and via a sneutrino (for the 
chargino) lead to di-lepton signatures of correlated flavor and charge with di-lepton masses significantly larger than 
for the X2 decay described above. While it seems very hard to disentangle the four corresponding kinematic end- 
points, the observation of heavy gaugino production is possible with a dedicated analysis in a part of the mSUGRA 
parameter space, if uiq is not too large. In particular, the highest mass end-point can be measured to a precision of 
approximately 4 GeV for the SPSla scenario with 100 fb" 1 of data. Its unambiguous identification needs additional 
information, however 



2.5. Measurement of the \'_! Spin 

The reconstruction of the supcrpartner spins is together with the determination of the their gauge quantum 
numbers the crucial test of SUSY. Spin information is carried by the scattering angle distribution of the primary 
pair of superpartners in their rest-frame. In hadron collisions, it is however hard to determine this frame due to 
the unknown longitudinal boost of the partonic initial state and due to the unobserved LSP's in the final state. It 
is however possible to exploit angular distributions of the supcrpartner decay products in decay chains to obtain 
spin information. A particular example for a measurement of the x° spin has been worked out, again exploiting the 
Ql — * q£ + £~Xi cascade (see Fig. 0J. Since the left squark always decays into a left-handed quark, the X2 becomes 
polarized. Since in its decay the xi emits a scalar right slepton, the corresponding so-called near lepton will carry 
the polarization information of the xi ■ Therefore one expects a charge asymmetry in invariant mass distribution of 
the quark and the near lepton, m(£ near q). For £ + q and t~ q, the tree-level differential form the m(£ near q) spectrum is 
dP/dm oc 4m 3 while for £~ q and £ + q it is dP/dm oc 4m(l — m 2 ). These tree level distributions are not experimentally 



observable since anti-quarks cannot be distinguished from quarks and the near lepton cannot be distinguished from 
the far lepton. If the original (anti)-squarks are (at least partially) produced by the qg — > qg and qg — > qg, the p.d.f. 
asymmetry due to the valence quarks in the protons can be exploited and more squarks than anti-squarks will be 
produced. The charge asymmetry from the far lepton is expected to be very small. Thus, forming both the invariant 
mass with the near and the far lepton dilutes the charge asymmetry but does not remove it. In Fig. [SJlcft) the 
reconstructed lepton-jet invariant mass distribution is shown for positive (squares) and negative (triangles) lepton 
charge. On the right, the resulting charge asymmetry is shown for 150 fb _1 |l4| . 
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Figure 8: Left: Distribution of the charged lepton - jet invariant mass distribution for positive and negative leptons after 
detector simulation. Right: Observable charge asymmetry in the charged-lepton mass spectrum after detector simulation for 
150 ftT 1 (from Q). 



3. SUSY AT THE LINEAR COLLIDER 



3.1. Experimental environment 



The International Linear Collider (ILC) is a projected electron positron collider at 500-1000 GeV center-of-mass 
energy with a luminosity of several 10 cm ^s" 1 H5lhdLJ- The energy will be tunable from the Z-pole up to the 
highest energy. Both beams can be polarized (e~: 90%, e + : 50-60%). The advantage of e + e~ collisions are the known, 
electro-weakly interacting initial state, the low level of instrumental and physics backgrounds and as a consequence 
the comparably low event rate which allows to record all collisions without any trigger requirements. Due to the 
high luminosity, unlike previous e + e~ machines, the ILC will have significant beam-beam interactions. These lead 
to a production of approximately 6xl0 10 low-energetic photons per bunch-crossing at 500 GeV. About 10% of the 
events will have center-of-mass energies below 95% of the nominal energy. Therefore the beam-strahlung spectrum 
has to be monitored continously with data (e.g. acollinearity of Bhabha events) and corrected for. Backgrounds from 
collisions of beam-strahlung photons have been studied and were found to be small except for the very forward region 
for which dedicated highly-granular calorimeters have to be build. 



3.2. Mass measurements 



At the ILC the masses of the color-neutral superpartners can be measured in two different ways. First in continuum 
production, kinematic end-points and energy spectra can be used to extract simultaneously the involved masses. 
Second, the measurement of the shape of the production cross-section for various processes near threshold allows for 
a very precise extraction of the sum of the produced superpartner masses. 

3.2.1. Sleptons 

Slcptons are pair-produced in the reactions 

e + e~ — > d'l' dj , P# &e I = e, /i, r and = L, R or 1,2] 

via s-channel 7 jZ exchange and t-channel x exchange for the first generation. 

As an example, in Fig. EJleft) the measurable energy spectrum of the muons from the process eje^, — ► jj^fi^ — * 
A t+ XiA t_ Xi i s shown [l8j. Events can be selected with negligible SM background. In particular background from 
W-pair production can be efficiently suppressed by choosing right-handed electrons in the initial state. SUSY 
backgrounds in this final state are generally small and can be suppressed in part by topological cuts. Due to the 
scalar nature of the smuons, the energy spectrum has a box shape. For the upper and lower end-points E+/_, the 
slepton and LSP masses can be determined as 



E_+E< 



Ve-Ej 




The situation is more complicated in the case of f slcptons due to the escaping neutrinos from r decay leading to 
a depletion of the of the upper end-point and to elimination of the lower endpoint. However, if the mass of Xi is 
known e.g. from smuon production, the shape and the upper end-point of the energy of the r decay products can 
still be used to extract a precise mass value for f\ as shown for the t — > 2>ttv t spectrum in Fig. fright). 




Figure 9: Left: Muon energy spectrum from the process e~£e R — > fifi — > ^Xif^ Xi Right: Hadron energy spectrum of the 
decay r — > 3wv T from the process eje^ — > fifi 

Alternatively, the slepton masses can be extracted from a threshold scan as shown in Fig. ^| for right selectron 
production both in e + e~ and e~e~ collisions and for right smuon production. With measurements at five ccntcr-of- 
mass energies with only 10 fb^ 1 per point a precision of 0(100 MeV) can be achieved. With this precision higher-order 
corrections and final width corrections have to be taken into account [nil ] . 

3.2.2. Charginos and Neutralinos 

Charginos and neutralinos are pair-produced 

e + e- -> xfxj = 1,2] (1) 

e + e- -> x°X° [i,i = l,...,4] (2) 
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Figure 10: cross-section for selectron and smuon pair production at threshold: &\e t 
eje^j — > fLnfiR (right). The error bars correspond to 10 fb _1 per point. 



ene R (left), e L e R — > e R e R (middle), 



via s-channel j/Z exchange and £-channel selectron or sneutrino exchange. The lightest chargino decays according to 
xf — > i^i/gXi either via an intermediate virtual or real boson or if kinematically possible via a real slepton. The 
second lightest neutralino decays according to £2 ~~ * ^ + ^~Xi either via a virtual or real Z boson or via a real slepton. 
In particular, if < m^o , invisible v decays may occur. For large mixing in the stau sector and for large values 
of tan (3 the f\ slepton is often much lighter than the other sleptons which can lead to a significant enhancement 
of r leptons in the chargino and neutralino final states. The production processes for f , X2 an d xt m ay therefore 
all lead to the same t + t~+ missing energy signature. Topological cuts and the use of polarized beams can help to 
disentangle the contributing SUSY processes. As in the case of sleptons, the chargino and neutralino masses can 
be measured from the lepton energy and mass spectra as well as from threshold scans. In the more difficult case of 
exclusive decays into r final states, a mass precision of a few GeV can be achieved in the continuum and 0.5 GeV 
from a threshold scan. Significantly better precision can be achieved if electron and muon final states are produced 
with sufficient rate [20j . 

3.2.3. Light Stop 

Although squarks are often too heavy to be produced at a 1 TeV LC, the light scalar top quark may be lighter 
than the other squarks and therefore accessible in the reaction e + e~ — > titi — > bx^bXi ~ * br + vx\br~ The 
final state consists of two 6-jets, two r's and missing energy. The energy spectrum of the &-jets can be used to 
reconstruct the stop mass provided the neutralino and chargino masses are known With a luminosity of 

1000 fb _1 the rate will be sufficient to achieve a mass resolution of 2 GeV. For a light scalar top quark, the decay 
chain e + e~ — > t\t\ > cxi ex? has also been studied. From a measurement of the production cross-section with 
opposite beam polarizations, a measurement of both mass and mixing angle can be inferred j^. 

3.2.4. Masses: Summary 

The achievable superpartner mass precision of the ILC for the SPSla scenario is summarized in Table [I] taken 
from 0. 



3.3. Quantum numbers, couplings and mixings 

Besides the precise measurement of the largest possible set of superpartner masses the measurement of quantum 
numbers, couplings, and mixings plays an important role in deciphering the supersymmetric model. In e + e~ collisions, 
due to the low background and the known initial state, various possibilities to extract quantum numbers and couplings 
exist. These range from the measurement of inclusive rates to the measurement of angular distributions in production 
and decay. 

3.3.1. Spin determination 

The fundaments of SUSY rely on the superpartners' spin differing by \ from their SM partners. It was shown in 
Sec. l2.5l that at the LHC a unique spin determination is quite involved. At the LC, the spins of the superpartners can 



Table I: Sparticle masses and their expected precisions in Linear Collider experiments, SPS la mSUGRA scenario (from [23j). 





m [GeV] 


Am [GeV] 


Comments 


xf 


176.4 


0.55 


simulation threshold scan , 100 fb _1 


xt 


378.2 


3 


estimate xtxf, spectra xt -> Zxt , Wx° 


X? 


96.1 


0.05 


combination of all methods 


X2 


176.8 


1.2 


simulation threshold scan X2X2, 100 fb _1 


X3 


358.8 


3-5 


spectra x° Zx?, 2 , X^Xa-XaX", 750 GeV, > 1000 fb" 1 


xi 


377.8 


3-5 


spectra xi - W X t, xlxl,X%xl, 750 GeV, > 1000 ib" 1 
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e~e~ threshold scan, 10 fb _1 
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e~e~ threshold scan 20 fb _1 
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simulation energy spectrum, 500 GeV, 500 fb _1 
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simulation energy spectrum, 400 GeV, 200 fb _1 
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estimate threshold scan, 100 fb _1 [24J 
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simulation energy spectra, 400 GeV, 200 fb _1 
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Figure 11: Distribution of the superpartner production angle for smuon (left) and charginos (right). The upper plots show 
both reconstructed solutions per event as points and the correct solution as full line; in the lower plots the combinatoric 
background is subtracted (from [rUn. 



be determined directly from the production angle distributions. The scalar leptons exhibit a sin 2 9 distribution which 
can be reconstructed up to a twofold ambiguity in smuon pair-production. After subtraction of the combinatorial 
background the angular distribution is very clean (Fig. lll| left) |l7j . The situation is more complicated for charginos 
and neutralinos which exhibit a forward-backward asymmetry in the production angle due to their mixed U(l) and 
SU(2) couplings and the additional t-channcl contribution. An example for chargino pair production is shown in 
Fig. 1111 right). The forward-backward asymmetry and in particular the left-right polarization asymmetry provide 
sensitive observables in order to disentangle the chargino and neutralino mixing matrices |^ . 

3.3.2. Chiral quantum numbers 

In SUSY, the chiral (anti-)fermions are associated in an unambiguous way to scalars, i.e. R <-> R and 
e \ R «-> eJj L . The four pair-production processes for left and right selectrons, e + e~ — > e^e^, e + e~ — > eje^, 



e + e~ — > e^e£, e + e _ — ► eje^ can be disentangled from their different dependence of the cross-section to polarized 
electron and positron beams. From Fig. 1121 it can be seen that e.g. e + e~ — > e^e^ and e + e _ — > eje^ have practically 
identical behavior of the cross-section as a function of the electron polarization but differ completely as a function of 
the positron polarization |26| . The t-channcl contribution to the production cross-sections is sensitive to the SUSY 
Yukawa coupling g(eex°) which is fundamentally related to the SM gauge couplings. The SU(2) and U(l) SUSY 
Yukawa couplin gs c an be determined to a precision of 0.7% and 0.2%, respectively with 500 fb _1 at 500 GeV in a 
SPSla scenario 1181. 



Selectron quantum numbers: unpolarised e+ Selectron quantum numbers: P(e-)=+90% 




P(e-) P(e+) 

Figure 12: Separation of the selectron pair e^e^ in e + e~ — > e\ R e~^ R may not be possible with electron polarization only 
(left); if, however, both beams are polarized, the RR configuration separates the pairs and the association of the selectrons to 
the chiral quantum numbers can be experimentally tested (from 1^3 1. 



3.3.3. Stau mixing angle 

In the third generation large mixing effects are expected. Again, beam polarization is vital to measure the stau 
mixing angle from measuring the cross-section for e + e~ — > f\T\ with polarized beams. This allows for a measurement 
of cos26f = — 0.84± 0.04 in an SPSla scenario for 500 fb _1 |2^. Further, independent information can be obtained 
from the measurement of the r lepton polarization in fx — > t\i — * n ^VXi- For known f mixing angle, this 
information can be used to gain sensitivity to tan (3. 

3.3.4. CP-violation in SUSY decays 

Many studies of SUSY at LHC and LC have focused on a (often constrained) MSSM with real parameters. The 
general MSSM Lagrangian however allows for CP-violating complex parameters. As an example, the U(l) gaugino 
mass parameter M\ and the Higgsino mixing parameter fi may have compl ex p hases eWi , <t>n- These phases may be 
accessed in angular correlations of the decay products in neutralino decay |28| . An experimental study of how well 
these correlations can be measured is still lacking. 

Similarly, CP-violating phases of the tri- linear couplings A T ,Ab,A t may be accessed in studying the branching 
fractions of the third generation sfermions an analyze them together with their masses and production cross sections 
in a global fit 0. 



3.4. Constraining Dark Matter 

The SUSY LSP provides an excellent candidate for dark matter. Recent measurements of temperature fluctuations 
of the cosmic microwave background by the WMAP satellite |3c| strongly constrain the SUSY LSP properties and 
therefore point to certain regions in the MSSM parameter space. Of particular interest for experimental studies 



at colliders is the co-annihilation region in which the neutralino annihilation is enhanced by the t-channel process 
\t — > T7 which contributes significantly only if the mass difference Am = m(r) — rn(xi) is small. The relic dark 
matter density depends critically on this mass difference. With the next generation of CMB experiments, in particular 
Planck, the DM density can be measured at the 2-3% level. It is therefore imperative to match this precision at 
colliders. 

If Am is small (typically below 10 GeV), the staus decay with small visible energy and the signature is only a 
few soft charged tracks accompanied by large missing energy. Two-photon background is becoming severe unless it 
can be efficiently vetoed by the detection of very forward scattered electrons. In the very forward region significant 
energy induced by beam-beam-interactions is deposited. This energy deposition in the most forward calorimeter, 
is shown in Fig. 1131 for two interaction region designs without (left) and with (right) a crossing-angle of the two 
incoming beams. The detection of energetic electrons is possible down to angles of 3.5 (5.7) mrad without (with) 
crossing angle if the calorimeter is finely segmented. 




The pair production of staus in the small-Am region has been studied in [3jJ and |22j for various MSSM parameter 
sets. With appropriate cuts, detection and a precise measurement of the f mass is possible down to Am ~3 GeV. The 
resulting precision on the prediction for the dark matter density ranges from 2 to 6%, depending on the f mass and 
on Am. This precision matches the anticipated precision of the Planck satellite of 2%. As an example the hadronic 
energy spectra for r decays from the process eje^ — > t\t\ — > t + x\ t ~Xi as shown after detector simulation and 
cuts together with the two-photon background for Am = 5 GeV (Model Point D' from [3). 



4. LHC/LC INTERPLAY 



The ultimate goal of measurements of the properties of superpartners at LHC and LC will be the extraction of the 
complete set of parameters of the low energy MSSM Lagrangian. The previous sections have already indicated the 
complementarity of the possibilities at the LHC (large mass reach for squarks and gluinos) and the LC (precision 
measurements of color-neutral part of spectrum) . In this section, we will discuss the additional benefit of simultaneous 
interpretation and possibly simultaneous data analysis if both machines run concurrently. This aspect has recently 
been studied in the international LHC/LC working group [2^]. Without aiming for completeness, a few examples of 
this interplay will be given in the following. 
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Figure 14: Hadron energy spectra E w of r — > 7r^ r , £" p of r — > and of r — > 37T!/ T decays from the reaction eJ j e R — > 
f i f i — » r + x? t~\i an d two-photon production assuming head-on collision. Model D', m,f l = 217.5 GeV, Am = 5.1 GeV, 
x/s = 600 GeV and £ = 300 ftr 1 



4.1. Joint analysis of superpartner masses 

It has been shown in Scction[2]that mass reconstruction of squarks and gluinos at the LHC suffers from the unknown 
masses of the lighter states, in particular the LSP and the sleptons. Under favourable circumstances, they can be 
reconstructed from a joint fit of various kinematic endpoints to moderate accuracy. However a strong correlation of 
e.g. the squark mass and the LSP neutralino mass remains (see Fig. ll5|) . A precise measurement of slcpton, chargino, 
and neutralino masses at the LC removes these correlations and improves the precision of squark and gluino masses 
considerably even in a situation where the latter are not directly accessible at the LC. The expected precisions on 
some of the masses in a SPSla scenario are shown Tab.[H] It should be noted that the precision on squark and gluino 
masses in the combined analysis is dominated by uncertainties on the hadronic energy scale of the LHC experiments, 
assumed to be 1%. 
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Figure 15: Correlation between left squark (left) and lighter sbottom (right) and the LSP mass at the LHC. 



Table II: The RMS values of the mass distribution in the case of the LHC alone, and together with the ILC measurement of 
m^o . Only the masses resulting from the edge analyses are given. All numbers in GeV. 



LHC LHC+ILC 



Am ; o 


4.8 


0.19 


Am; 


4.8 


0.34 


Am ; o 


4.7 


0.24 


Am,g L 


8.7 


4.9 


Am h 


13.2 


10.5 



4.2. ILC mass predictions for LHC searches 

At the ILC, the SUSY parameters which govern the chargino-neutralino sector, i.e. the U(l) and SU(2) gaugino 
mass parameters Mi,M2, the Higgsino mixing parameter /i and tan/3 can be uniquely and precisely extracted from 
the measurements of masses and polarized cross-sections of the lightest gauginos, Xi)X2i an d Xi ■ These parameters 
can in turn be used to predict the masses of the heavier members x!j;X4 an d X2 within the framework of the 
general MSSM. With such a mass prediction, the search for the (often small) signals of the heavy gauginos can be 
substantially facilitated. The search for an edge in the di-lepton mass spectrum becomes transformed into a single 
hypothesis test, thus inceasing the statistical power of the LHC data for this specific hypothesis. 

Should in turn the predicted state be observed at the LHC, its measured mass can be fed back into the SUSY 
parameter analysis and considerably improve the achieavblc precision. 

Furthermore, the comparison of the predicted xl mass from LC measurements and the measured mass at the LHC 
allows to test for models beyond the MSSM. This has recently been shown J3] for a NMSSM, where in contrast to 
the MSSM five neutralinos are predicted with a mass spectrum not satisfying the MSSM relations. 



4.3. Global MSSM fits and reconstruction of the theory 

Ultimately, the goal is to measure the complete set of electro-weak scale parameters of the SUSY Lagrangian. 
At tree level, the parameter determination can proceed sector by sector, e.g. the chargino sector is completely 
determined by the three parameters Af 2 , /i, tan/?. However, with the anticipated precision of ILC measurements, 
higher order corrections to masses, cross-sections, and branching ratios are not negligible. At loop-level in principle 
every observable depends on the full set of SUSY parameters. An analytic procedure to extract the Lagrangian 
parameters from data is no longer possible. Instead, a global fit of the Lagrangian parameters to the complete set of 
SUSY observabales at LHC and LC will be necessary. 

Two programs, SFITTER 35] and Fittino [3| have been developed to achieve this goal. As an example, a 
recent result from Fittino is explained here. A MSSM with 19 free parameters has been chosen as the theoretical 
basis. It is derived from the full MSSM but assuming real couplings, flavour diagonal sfermion mass matrices and 
universality of the soft SUSY breaking parameters of the first t wo g enerations. For the theoretical predictions for 
the observables as a function of the parameters the SPHENO [37j program has been used. This code includes 
higher-order corrections whereever they have been calculated. The program is interfaced via the 'SUSY Lcs Houches 
Accord' (SLHA) [ssj . a convention to exchange SUSY parameter information between various programs in a coherent 
way. Therefore it is possible to replace the actual SUSY code behind Fittino in order to perform comparisons. As 
simulated measurements, the masses measureable at LHC and LC as explained in the previous sections have been 
used. In addition polarized cross-section measurements at 500 and 1000 GeV LC have been input with accuracies 
based on estimates from the expected statistical errors, but including systematic errors of 1% in each case. For the 
true values of the parameters, the SPSla scenario has been chosen but nowhere in the fit any assumptions on the 
particular SUSY breaking mechanism have been made. Special care was given to a fitting strategy which does not 



Table III: The Fittino SPSla fit result. The left column shows the assumed SPSla values, the middle column represents the 
result of the intermediate fit without m top , and the right column shows the result of the final fit. All SPSla input values of 
the parameters are reconstructed. 



Parameter 


SPSla value 


Tree-level 
estimate 


Final fit result 


tan 3 


10.0 


9.97 


10.0 ± 0.3 


u 


358.64 GeV 


354.4 GeV 


358.6 ± 1.1 GeV 


X T 


-3837.23 GeV 


-3533.0 GeV 


-3837.2 ± 131.0 GeV 


M e „ 


135.76 GeV 


150.2 GeV 


135.76 ± 0.39 GeV 


T R 


133.33 GeV 


141.0 GeV 


133.33 ± 0.75 GeV 


Me, 


195.21 GeV 


202.7 GeV 


195.21 ± 0.18 GeV 


Mfr 


194.39 GeV 


206.6 GeV 


194.4 ± 1.18 GeV 




-506.388 GeV 


-43.5 GeV 


-506.4 ± 29.5 GeV 


V 

bottom 


-4441.0 GeV 


-3533.0 GeV 


-4441.1 ± 1765 GeV 


Mr 
»B 


528.14 GeV 


567.3 GeV 


528.2 ± 17.6 GeV 


Mr 
b R. 


524.718 GeV 


566.0 GeV 


524.7 ± 7.7 GeV 




530.253 GeV 


566.9 GeV 


530.2 ± 19.1 GeV 




424.382 GeV 


373.7 GeV 


424.4 ± 8.54 GeV 


Mu L 


548.705 GeV 


581.3 GeV 


548.7 ±5.2 GeV 


M- H 


499.972 GeV 


575.4 GeV 


500.0 ±8.1 GeV 


Mi 


101.809 GeV 


99.07 GeV 


101.81 ±0.06 GeV 


M 2 


191.7556 GeV 


195.08 GeV 


191.76 ±0.10 GeV 


M 3 


588.797 GeV 


630.5 GeV 


588.8 ± 7.9 GeV 


m A 


399.767 GeV 


399.8 GeV 


399.8 ± 0.71 GeV 


m top 


174.3 GeV 


174.3 GeV 


174.3 ± 0.34 GeV 



make use of any a-priori knowledge of the parameters. In particular, suitable start values for the parameters have 
been estimated from the measurements using tree-level relations between various observable and parameters sector by 
sector. In order to yield a correctly converging fit an iterative procedure has to be applied which does not leave free 
all parameters from the beginning. Only after this preparatory phase the full fit can be performed with all parameters 
left free. The example result for SPSla is shown in Table ITTTl The importance of higher-order corrections can be seen 
from comparing the tree-level parameter estimates with their final values. It has been verified that the obtained fit 
errors are consistent with the fluctation of the fitted parameters in many repeated measurements. It should be noted 
that neither LHC nor LC input alone can constrain the assumed model enough to yield a converging fit. It was also 
shown, that due to the mutual influence of many parameters in determining a single observable, wrong assumptions 
of un-fitted parameters lead to wrong central values for the fit parameters. The lesson from this excercise are: first, it 
is possible, at least for a favourable scenario like SPSla, to extract the complete electro- weak scale Lagrangian from 
the future measurements of LHC and LC without strong assumptions on the particular SUSY breaking scenario. 
Second, the fact that results from both LHC and LC are necessary to obtain this result nicely shows the strong 
interplay between both machines. Third, the achievable experimental precision clearly requires the knowledge of 
theoretical predictions beyond leading order. The definition of a clear scheme to extract well-defined par ameters at 
higher orders is currently being worked out in the Supersymmetry Parameter Analysis (SPA) project |39j. 

The extracted parameters of the electro-weak scale MSSM Lagrangian can then be extrapolated to high (GUT, 
Planck) scales (Fig. I16l in order to determine distinct patterns of unification and reconstruct the underlying funda- 
mental theory of SUSY breaking |4cj . 
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Figure 16: Renormalization group evolution of the MSSM parameters determined from LHC and LC measurements to the 
GUT scale. 
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